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Noise and instability of an optical lattice clock 
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We present an analysis of the different types of noise from the detection and interrogation laser in 
our strontium lattice clock. We develop a noise model showing that in our setup quantum projection 
noise-limited detection is possible if more than 130 atoms are interrogated. Adding information 
about the noise spectrum of our clock laser with sub-10~^® fractional frequency instability allows 
to infer the clock stability for different modes of operation. Excellent agreement with experimental 
observations for the instability of the difference between two interleaved stabilizations is found. We 
infer a clock instability of 1.6 x 10“^®/-y/r/s as a function of averaging time r for normal clock 
operation. 

PACS numbers: 42.62.Eh, 32.30.-r, 37.10.Jk 


I. INTRODUCTION 

Optical clocks in general and lattice clocks in particular 
have shown outstanding fractional frequency instabilities 
of about 2x 10“^®/y^T/s [IHS], averaging down to the low 
10“^®-regime. This development has been made possible 
mainly by significant improvements on the frequency sta¬ 
bility of the interrogation lasers [H]. The improvement 
of clock stability has a strong impact on the determina¬ 
tion of the clock’s accuracy since systematic frequency 
shifts can be evaluated with higher accuracy in reason¬ 
able time. Furthermore, the stability of a clock also de¬ 
termines how practical it is for actual measurements and 
if, e.g., temporal variations of signals [7] can be observed. 

The instability of a clock stems from different sources 
of noise contributing to the error signal detected in the 
clock cycle, i.e. the estimated frequency offset of the clock 
laser from the atomic transition, and from the Dick effect 
[8] due to the non-continuous interrogation of the atoms 
by the clock laser [5]. The Dick effect originates mostly 
from the unobserved clock laser fluctuations during the 
dead-time of the clock cycle. 

In high-performance clocks the reference transition 
is coherently interrogated by either Rabi or Ramsey 
schemes on a frequency-sensitive slope of the spectro¬ 
scopic signal. In the detection, the quantum superpo¬ 
sition state of each individual atom is projected to the 
ground or excited state, leading to the fundamental noise 
limit given by the quantum projection noise (QPN) |I0j . 
An increase of atom number N will reduce the influence 
of the projection noise as it scales with \/N while the sig¬ 
nal scales with N. However, in practice the atom number 
is limited by other factors as, e.g., collision shifts [IIHIll- 
Ideally, the QPN should be the dominant noise contribu¬ 
tion in the measurement of the frequency offset of the 
local oscillator from the atomic line. 

Squeezed states m and entanglement of the atoms 
m can overcome the projection noise limit. It must, 
however, be noted that at present even the most sta¬ 
ble optical clocks are still limited by the aliasing of laser 
noise, i.e., via the Dick effect. 

With this background it becomes obvious that a de¬ 


tailed understanding of the noise sources present in the 
experiment is not only essential for optimizing the clock 
stability and interrogation strategy, but also for judging 
the necessity to implement squeezing and entanglement 
methods to improve the clock. 

We therefore give a rigorous analysis of the detection 
noise contributions in our Sr lattice clock [n! (Sec.|n|. 
The noise model is further supplemented by a contri¬ 
bution of the clock laser noise [5] via the Dick effect 
(Sec. III). Combining both, we show that the clock insta¬ 
bility observed during evaluations of systematic effects in 
interleaved stabilization mode is well reproduced and a 
clock instability of 1.6 x ip-i s/ ^/ t/s can be inferred for 
regular clock operation (Sec. IV), which is governed by 
the excellent stability of our clock laser [6] and compet¬ 
itive with the best values of 2.2 x 10“^®/-\/t/s [T] and 
3.2 X 10“^®/-\/r/s [2] published to date. 


II. DETECTION NOISE OF THE LATTICE 
CLOCK 

In this section, the individual noise contributions af¬ 
fecting the measured excitation probability are discussed. 
In the final part of the section, their influence on the clock 
signal is modeled and compared with experimental obser¬ 
vations. For the sake of simplicity, we express all noise 
amplitudes in arbitrary units of the data acquisition sys¬ 
tem labeled ‘counts’. 

Figure shows a simplified level scheme of the stron¬ 
tium atom. Details for laser cooling and trapping of 
strontium atoms have been described in previous pub¬ 
lications [TM9]. To derive the spectroscopic signal of 
the lattice clock after the interrogation of the ®^Sr atoms 
by the clock laser, atoms in the ground state (^^o) are 
excited with a resonant laser beam in standing-wave 
configuration on the strong 461 nm cooling transition, 
— ^Pi- The fluorescence is observed by a photomul¬ 
tiplier tube in current mode. Its signal is amplified and 
digitized with an analog-to-digital converter of the data 
acquisition computer. This signal g is, apart from an 
offset o, proportional to the ^5'o ground state atom num- 
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FIG. 1: Simplified level scheme of the strontium atom. 
Arrows indicate the transitions used for cooling and spec¬ 
troscopy and their associated wavelengths. 


to a steeper slope of 
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A. Electronic noise 

The electronic noise was measured by running the stan¬ 
dard detection sequence without atoms or laser light be¬ 
ing present. The observed noise is a sum of amplifier, dig¬ 
itizer, and dark current noise. Such noise will be present 
in the three signals, g, e, and o, but not in d. The 
standard deviation of the signal is deiec = 0.82 counts, 
while the amplitude of the electronic offset signal S'eiec ~ 
18 counts. 


ber. The radiation pressure removes the atoms from the 
detection volume within 20 ms. 

Similarly, the signal e from atoms in the excited clock 
state {^Pq) is detected after the atoms have been op¬ 
tically pumped within few 100 /iS to the ^Pi state by 
two lasers resonant with the ^Pq — ^Si and ^P 2 — ^S'l 
transitions, from where they decay further to the ground 
state. In addition, the offset o in both g and e due to 
stray light, multiplier dark current, and electronic offset 
is measured via a third and final detection pulse after 
removing any remaining atoms. Therefore, each of the 
signals 5 , e, and o is a sum of different contributions 
from fluorescence (S’auo)) stray light (S'stray)j and elec¬ 
tronic background (S'eiec)- From these three signals, the 
atomic excitation probability Pe is estimated by 


This normalization removes substantial noise on g and e 
due to shot-to-shot fluctuations of the atom number N. 
Here, we replaced the measured offset o by its running 
average b, since the variation of the actual background is 
small on timescales of several interrogation cycles; thus 
the additional noise of Pe due to the shot-to-shot noise of 
the offset measurement is suppressed and does not need 
to be considered in the noise model. 

The connection between the excitation probability, in 
particular its noise, and the corresponding frequency ex¬ 
cursion of the interrogation laser is given by the slope 
of the spectroscopic signal. For the case of probing the 
atomic resonance line at a half-width point [S] using Rabi 
interrogation with a pulse of length the slope is 

^ « ±27r • 0.30 • (2) 

av 

with its sign depending on which side of the resonance 
is probed. A Ramsey interrogation scheme with a free 
precession time TRamsey and short excitation pulses leads 


B. Photon shot noise 

During the detection of the atomic fluorescence, a fi¬ 
nite number of photons is collected by the photomulti¬ 
plier tube. Thus, the signals g and e will suffer from a 
shot noise contribution. In order to quantify this contri¬ 
bution, we have investigated the photon shot noise using 
a flashlight as a shot noise-limited light source. For the 
observed white phase noise, the first point of the Allan 
deviation of the recorded data is equal to their standard 
deviation; we use this value as a measure of the noise to 
remove the influence of slow intensity variations of the 
flashlight. The measurements were performed for differ¬ 
ent signal amplitudes. The electronic background S'eiec 
(discussed in the previous subsection) must be removed 
to extract the actual ‘fluorescence’ signal Sfluo, be. the 
signal stemming only from detected photons. The re- 
sultiM ‘fluorescence’ noise contribution Csn is shown in 
Fig.[^ The expected VSfluo-dependence of the noise am¬ 
plitude (Tsn is well reproduced. From a fit we find 

^s„(Sfluo) = 0.59(2)v^, (4) 

which corresponds to Sfluo = 0.35(3) counts per detected 
photon. Fluorescence shot noise will be present in the 
signals g and e with an amplitude depending on the flu¬ 
orescence contribution Sfluo to these signals, as e.g. the 
electronic offset Seiec is not subject to shot noise. 


C. Detection laser intensity noise 

Intensity fluctuations of the detection laser at the posi¬ 
tion of the atoms will show up on the fluorescence signal 
as long as the transition is not strongly saturated. We 
avoid high intensity, since the radiation pressure-induced 
heating of the atoms reduces the reasonable interaction 
time and thus the detected signal. This loss of signal is 
not compensated by the higher photon scattering rate, 
and the detection would be less efficient. 
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FIG. 2: Measured fluorescence shot noise levels Usn (full cir¬ 
cles) as a function of the detected signal 5fiuo along with a fit 
according to Eq. (solid line). 


Intensity fluctuations may arise from both power and 
pointing instability of the detection laser. We measured 
the power fluctuations of the laser in a similar procedure 
as for the shot noise measurements, except that we used 
stray light from the detection laser instead of a flashlight. 
From these measurements we have found data that is 
very similar to the one presented in Fig.[^ In particular, 
we have observed no significant noise contribution with 
a linear dependence on laser power and thus conclude 
that laser power noise on short time scales, which would 
result in such a contribution, can be neglected. This is 
corroborated further by direct measurements of the laser 
power. The optical setup of the detection beam, which is 
delivered by a fiber and collimated to a diameter of about 
2 mm, comprises only a short free-space path on the order 
of 50 cm; we have analyzed its pointing instability and 
found it to be negligible on relevant time scales. We also 
note that the shot noise of the detection beam with a 
power of about 400 /rW is negligible. 

Laser stray light (Astray) adds a shot-noise contribution 
to the signals g and e according to Eq. whereas the 
similar noise contribution to the offset is suppressed by 
the use of o. 


D. Detection laser frequency noise 

The 461 nm detection laser beam is derived from a 
frequency-doubled diode laser system that also produces 
the laser beams for laser cooling and Zeeman slowing. 
The fundamental frequency of the laser system is stabi¬ 
lized to a high-finesse, 10 cm long ultra-low expansion 
glass (ULE) resonator. With its resonance line width of 
less than 100 kHz, we estimate in-lock frequency fluctua¬ 
tions on the kilohertz level from the in-loop error signal. 
Compared to the 32 MHz line width of the de¬ 

tection transition, frequency noise of the detection laser 


FIG. 3: Excitation probability noise atpe = 0.5 multiplied 
by the signal amplitude g -I- e — 2o as a function thereof (full 
circles). Atoms are prepared in a superposition of and 
^Po states by a frequency-insensitive resonant 7r/2-pulse. The 
dashed line shows a fit of Eq.|^with a factor of proportionality 
close to the experimentally expected one (see text). 


does not contribute significantly to the detection noise. 


E. Quantum projection noise 


As mentioned in Sec. [Tj QPN is the fundamental noise 
that ideally should dominate all other noise contributions 
in the detection. It depends on the atom number N and 
the excitation probability Pe [ID] as 


O'QPN(Pe) = VPei^ - Pe)N. (5) 


The signal-to-noise ratio Af/crqpN can thus be improved 
by increasing the atom number. N is proportional to 
g+e—2o, where the constant of proportionality to convert 
from ‘counts’ to atom number is about one atom per 
count. This factor depends on the actual alignment of 
the experiment, the frequency and power of the detection 
laser, and other parameters. 


We measured the noise ap^ of the excitation proba¬ 
bility of samples of atoms with pe = 0.5 and different 
atom numbers N to investigate whether we can achieve 
QPN-limited detection in our experiment. In order to 
become insusceptible to frequency noise of the interroga¬ 
tion laser, we prepared a coherent superposition state by 
using a resonant 7r/2-Rabi pulse instead of Rabi excita¬ 
tion with a TT-pulse at a half-width frequency detuning 
as usually applied in a stabilization sequence. A pulse 
length of T ^/2 = 10.5 ms was chosen. The observed noise 
CTp^ is plotted in Fig.[^ To which extent quantum projec¬ 
tion noise is dominant in our setup and if residual laser 
noise affected the data will only become apparent with 


the combined noise analysis in Sec. HE 
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FIG. 4: Experimental noise of the excitation probability with 
suppressed sensitivity to laser frequency noise (7r/2-pulses, 
full circles) and estimated individual noise contributions. The 
QPN is based on a detection efficiency of 0.65 atoms/count. 
The shot noise calculation includes the noise of a typical back¬ 
ground of 70 counts due to detection laser stray light. The 
green curve shows the summed noise according to Eq. The 
good agreement with the experimental data demonstrates the 
completeness of the model and the absence of laser noise in 
the data. 


F. Detection noise model 

Having quantified the individual noise sources, we now 
develop a noise model to combine them and verify if the 
observed detection noise (Fig. is fully described and 
whether QPN is the dominant noise source. 

The noise model assumes independent contributions to 
tjp^ from QPN and other sources of noise. QPN is han¬ 
dled separately because it leads to anti-correlated noise 
in g and e. For the other noise sources, independent con¬ 
tributions from g and e, but not from b, are considered. 
For these contributions, the model uses the derivatives of 
Pe (Eq. 0 with respect to g and e. Thus the total noise 
is given by 



( 6 ) 


with the derivatives 

dpe d — e 

dff {g + e- 2b)2 

dpe _ g-o 

de {g + e — 2o)^ 


FIG. 5: Sensitivity function w{t) of interleaved stabilizations, 
shown for the cases of lumped (solid line) and distributed 
(dashed line) arrangements, as described in the text. Vertical 
grey lines indicate multiples of Tc. The typical parameter 
values given in Sec. |III| were used. For a single stabilization, 
w{t) consists only of the first four cycles shown for the lumped 
interleaved arrangement. 


in reasonable agreement with an independent calibration 
derived from absorption measurements. 

From this analysis we conclude that detection is limited 
by QPN for more than 200 counts, or 130 atoms. For 
300 atoms this would lead to a frequ ency stability of our 
lattice clock of less than 6 x 10 “^’^/a/t/s, where we made 
use of Eqs. and and used a realistic interrogation 
time T,r = 640 ms with a cycle time of T(. = 1 s [T71 HD] . 
However, this estimate does not take into account the 
degradation of the stability due to the Dick effect, i.e. 
the aliasing of high frequency laser noise due to the non- 
continuous interrogation of the atoms. 


III. LASER NOISE AND DICK EFFECT 


After establishing the noise limitations on the optical 
clock instability from the atomic side, we will now show 
what frequency instability is introduced by the interroga¬ 
tion laser. Because of the non-continuous interrogation of 
the atoms, a stability degradation due to aliasing of noise 
- the so-called Dick effect [5] - is expected. Its contribu¬ 
tion, CTp,Dick('r), to the fractional clock instability can be 
calculated by 


11 

^y.Dick('^) = H '^!/(^/'^) IS/cT 

^ l50l 


( 8 ) 


In Fig. we summarize the contributions to the de¬ 
tected excitation probability noise as a function of the 
total atom number as expressed by g + e — 2d. We see 
that the total noise agrees with the observations from 
Fig. i very well, when we use a conversion factor of 
0.65 atoms/count to calculate the QPN in Eq. This is 


where T = nTc, n is the number of interrogations in each 
stabilization cycle. The cycle duration Tc consists of the 
preparation time Td and the interrogation time T,^- 
In Eq.|^ the laser’s single-sided power spectral density 
Sy is evaluated at multiples of the inverse of the total du¬ 
ration T of a complete stabilization cycle weighted with 
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FIG. 6 : (a) Spectral power density of frequency fluctuations Si, = v^Sy (solid line) and Fourier components | 5 fe/Pnorml (^rcles) 


for k > 0 for interleaved stabilizations for parameter values as given in Sec. Ill (b) Partial sums o f o'^ .pick)'’") (Eq. as a 


function of the Fourier frequency /max = fcmax/F at r = 1 s and for the parameter values given in Sec. Ill Shown are the cases 
of interleaved stabilizations using either a lumped (solid circles) or a more stable distributed (open circles) probing sequence, 
as described in the text, and single stabilization using Rabi (solid squares) and alternatively Ramsey (solid triangles, same 
interrogation time with 7 r/ 2 -pulses of 100 ms length) interrogation. 


coefficients gk, which are the Fourier components of a sen¬ 
sitivity function w(t)^ and normalized by the dc Fourier 
coefficient g^. The sensitivity function w{t) describes the 
response of the detected excitation probability to inter¬ 
rogation laser frequency changes [5] 

5pe = X / 2tt Sv (t) ■ w (t) dt. (9) 

2 Jo 

The dc Fourier coefficient g^ = 1 /tt dpe/dv in Eq. [pre¬ 
lates the change in excitation probability to a constant 
change in laser frequency. 

For a TT-pulse of duration centered at t = TtiI2 and 
a detuning A « iO.dO/Tn- of the interrogation laser, i.e., 
to the half-maximum point of the resonance, it can be 
found (equation 11 in [5]) that 

{ sin^ J cos t) 

X [(1 — cos 1 / 2 ) sin 

- 1-(1 — coslli) sinn 2 ] during pulse, 

0 elsewhere 

with 

■d = ^ — arctan(2T^A) (10) 

and 

fli = ^'/iT(2T^x ^ 

^ TT 

fl2 = ^yi + (2T,A)" X 


During clock operation, each stabilization cycle of the 
clock consists of four distinct interrogation (and prepara¬ 
tion) sequences that interrogate the mp = ±9/2 compo¬ 
nents on both slopes m- Thus, a clock laser frequency 
correction is applied every 4Tc, and w{t) is a fourfold 
series of w single (0 (Fig-[p)- 

We use an interleaved stabilization scheme to evaluate 
most systematic shifts. The parameters of interest are 
alternated after each full clock stabilization cycle; sep¬ 
arate digital servos correct the clock laser frequency in¬ 
dependently for each configuration and generate a differ¬ 
ence frequency signal. As the laser frequency noise is in 
general correlated between the interleaved interrogations, 
the Dick effect cannot be simply treated independently 
for each of the stabilizations. In analogy to the Dick ef¬ 
fect for clock operation, a Dick effect can also be derived 
for the difference between the two interleaved stabiliza¬ 
tions. The sensitivity function is composed of eight rep¬ 
etitions of Wsingie(i)) where the sign is reversed between 
the first and last four (Fig.jp); a possible improvement to 
this lumped arrangement is discussed in Sec. [^ below. 
The conversion coefficient between the signal, i.e., the 
difference of excitation probabilities, and laser frequency 
changes is no longer given by go^ but by the response 
to a unit frequency difference between both stabilization 
settings, which is 

5 norm=y ^\w{t)\dt. ( 11 ) 

Thus, the instability of the difference between the two 
interleaved stabilizations can be calculated from Eq. jp 
with the above sensitivity function and replacing g^ by 

^norm* 























6 



10 ° 10 ' 10 ° 10 ° 10 ° 
averaging time x (s) 

FIG. 7: Allan deviation of the difference between two in¬ 
terleaved stabilizations with = 0.64 s, Td = 0.54 s and 
N ~ 360 atoms (full circles). The fractional instabilities (Jy(r) 
inferred from our analysis for interleaved stabilizations (solid 
blue line, 4.7 x 10“^®/-y/r/s) and pure clock operation (dashed 
red line, 1.6 x 10“^®/'\/t/s) are shown for comparison. 


To calculate from Eq. the laser noise 

spectrum is required. We have measured Sy via three- 
cornered-hat comparisons with other lasers [S] . The noise 
spectrum is shown together with the gk for fc > 0 in 
Fig. |6]t. 

Having now all ingredients, the sum in Eq. can 
be evaluated for r = 1 s and typical parameter values 
(Tn- = 0.64 s, Td = 0.54 s, and Tc = T,r + Td = 1.18 s). 
To visualize the contributions to cxy^Dick from different 


Fourier frequencies, the partial sum of Eq. 
frequency /max = ^max/T is plotted in Fig. 
that for a single stabilization, a. 


y,Dick V 


up to the 
We see, 
is dominated by 
gi, whereas higher frequencies between 1 Hz to 10 Hz 
contribute significantly for the interleaved stabilization. 
In fact, because of the symmetry of the single stabiliza¬ 
tion, all gk where k is not a multiple of 4 are zero. Thus, 
the Dick effect is exactly the same as for a single inter¬ 
rogation cycle of length Tc. 


IV. LATTICE CLOCK INSTABILITY 

To estimate the instability of our lattice clock, the re¬ 
sults from Secs. |H and |lTT| were combined; the instabil¬ 
ities due to detection noise and the Dick effect have to 
be added in quadrature. For the difference of two inter¬ 
leaved stabilizations, the detection noise contribution is 
\/2 times that of a single stabilization discussed in Sec.[ll| 
due to the two independent stabilizations, whereas the 
Dick effect for the difference is accounted for by the com¬ 
bined sensitivity function discussed in Sec. |HI[ 

From our analysis we infer a combined instabil¬ 
ity of 4.7 X lo^is/y'r/s, with contributions of 4.6 x 
10 “^®/a/t/s from the Dick effect and 6 x 10“^’^/y^T/s 


for each of the interleaved stabilizations from detection 
noise. Experimentally, an instability very close to that 
inferred from our analysis is found. Figure shows the 
Allan deviation of the difference signal observed in an 
interleaved stabilization for the experimental parameters 
given in Sec. |HI| and N « 360 atoms. Thus, we conclude 
that the combined analysis of detection and laser noise 
constitutes a very good description of our optical lattice 
clock. 

The stability of interleaved stabilizations can be im¬ 
proved by modifying the sequence to use a distributed 
instead of a lumped arrangement, i.e. interleaving the in¬ 
terrogation sequences of each stabilization as well rather 
than arranging them in blocks of four, as shown in Fig.l^ 
We estimate that a combined instability of 3.4 x 10“^ 
in 1 s can be achieved. 

Based on the previous hndings, we can infer the insta¬ 
bility in clock operation with typical atom numbers of 
N « 300 atoms. Experimentally, we cannot determine 
this quantity directly, since we are lacking sufficiently 
stable oscillators to compar e to . With this analysis we 
find cry{T) = 1.6 x 10“^®/-\/t/s dominated by the Dick 
effect-induced instability of ay Dick(T) = 1.5x10“^®a/t/s 
(Fig.0. A slight reduction of instability can be achieved 
by Ramsey interrogation instead of Rabi pulses (Fig. |^. 


V. CONCLUSION 

We have given a full analysis of the detection noise of 
our Sr lattice clock. In combination with a Dick effect 
analysis of our new clock laser [5], we are able to give a 
full description of our clock’s observed frequency insta¬ 
bility when differentially evaluating systematic frequency 
shifts. Since the agreement between model and observa¬ 
tion is excellent, we are convinced that the instability of 
CTy(r) = 1.6 X 10“1®/.^t/s we determine for normal clock 
operation is realistic. This is an exceptionally small in¬ 
stability, better than the so far published instabilities of 
optical clocks m- 

Although the residual noise of our clock laser [5] has 
been crucial for achieving this result, we also conclude 
that, even with one of the most advanced interrogation 
lasers available today and an efficient preparation scheme 
with a duty cycle of more than 50 %, the clock instability 
is, already at very small atom numbers, limited by the 
Dick effect and thus the clock laser. This means that at 
present advanced squeezing or entanglement methods are 
not yet worthwhile. 

First, the Dick effect-induced instability must be ad¬ 
dressed. Ramsey interrogation offers a favorable sensitiv¬ 
ity function as compared to a Rabi scheme, especially for 
high duty cycles [S]. Moreover, reducing the clock laser 
frequency noise between 1 Hz and 10 Hz, which is caused 
by seismic perturbations of the reference cavity, e.g. by 
active feed-forward [2T], could bring the clock instability 
to below 10“^® in one second. Finally, the duty cycle of 
clock laser interrogation could be increased: Nondestruc- 
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tive detection methods [22], which use, e.g., the phase- 
shift imprinted onto an off-resonant detection beam, al¬ 
low keeping the majority of cold atoms in the lattice and 
thus reduce the required preparation time. Furthermore, 
few independent physics packages can use the same clock 
laser for interrogation to achieve a dead time-free ob¬ 
servation of the laser frequency, thereby eliminating the 
Dick effect altogether. The detection noise presented in 
Sec. 1^ would also be reduced by a factor of owing to 
the dual interrogation per cycle. Such a setup leads to an 
ultimate instability of 4 x 10 “^^/yjrjs for the typical pa¬ 
rameters discussed here and allows quantum projection 
noise-limited detection. 
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